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Abstract. We give a comprehensive study on regularized approximate 
electromagnetic cloaking via the transformation optics approach. The 
following aspects are investigated: (i) near-invisibility cloaking of passive 
media as well as active/radiating sources; (ii) the existence of cloak- 
busting inclusions without lossy medium lining; (iii) overcoming the 
cloaking-busts by employing a lossy layer outside the cloaked region; 
(iv) the frequency dependence of the cloaking performances. We address 
these issues and connect the obtained asymptotic results to singular 
ideal cloaking. Numerical verifications and demonstrations are provided 
to show the sharpness of our analytical study. 



1. Introduction 

A region is said to be cloaked if its contents together with the cloak 
are invisible to certain noninvasive wave detections. Blueprints for making 
objects invisible to electromagnetic waves were proposed by Pendry et al. 
[18] and Leonhardt [11] in 2006. In the case of electrostatics, the same idea 
was discussed by Greenleaf et al. [8] in 2003. The key ingredient for those 
constructions is that optical parameters have transformation properties and 
could be push- forwarded to form new material parameters. The obtained 
materials/media are called transformation media, which we shall further 
examine in the current work for cloaking of the full system of Maxwell's 
equations. 

The transformation-optics-approach-based scheme proposed in [8, 18] is 
rather singular. This poses much challenge to both theoretical analysis 
and practical fabrication. In order to avoid the singular structures, several 
regularized approximate cloaking schemes are proposed in [6, 9, 10, 12, 22]. 
The works [6] and [22] are based on truncation, whereas in [9, 10, 12], the 
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'blow-up-a-point' transformation in [8, 18] is regularized to be the 'blow- 
up-a-small-ball' transformation. The performances of both regularization 
schemes have been assessed for cloaking of acoustic waves to give successful 
near-invisibility effects. Particularly, in [9], the authors show that in order 
to 'nearly-cloak' an arbitrary content, it is necessary to employ an absorbing 
('lossy') layer lining right outside the cloaked region. Since otherwise, there 
exist cloaking-busting inclusions which defy any attempts of cloaking. This 
idea of introducing a lossy layer has recently been intensively investigated for 
approximate acoustic cloaking (see [15, 16]), whose behaviors are now much 
well understood. However, very little progress has been made in the study 
of approximate EM cloaking for full Maxwell's equations due to the more 
complicated structure of Maxwell's equations. This is the main concern of 
the present article. 

We have considered both the 'truncation-scheme' and the 'blow-up-a- 
small-ball-scheme' for approximate EM cloaking. However, in our study, the 
two regularization schemes have the same performances for near-invisibility, 
so we focus our exposition on the latter one. Based on a model problem, the 
following aspects on the approximate EM cloaking are addressed in detail. 

(i) The near-cloak of EM waves for both passive media and active / radiating 
sources. For approximate cloaking of passive media, the near-cloak is shown 
to be within p 3 of the singular ideal cloaking, where p is the regularization 
parameter. Whereas if there is a delta point source present in the cloaked 
region, the near-cloak is shown to be within p 2 of the perfect cloaking. That 
is, we could still achieve the near-invisibility effect, but with one order reduc- 
tion on the approximation. Compared to the near-invisibility assessments 
in [6, 9, 12, 15, 16, 22] for approximate acoustic cloaking (which is of 0{p) 
when spatial dimension is 3; and 0(1/ In p) when spatial dimension is 2), 
the performances for near-cloak of EM waves are much better. We point out 
that the study in [9, 12, 15, 16, 22] lacks the analysis on the approximate 
cloaking when there is an active/radiating source present inside the cloaked 
region. Another rather interesting observation we would like to make is 
that in [5], it is shown one cannot perfectly cloak an H~ l -source inside the 
cloaked region since otherwise there would be a conflict with certain 'hidden' 
boundary conditions of the finite energy solutions underlying the singular 
ideal EM cloaking, but our analysis here shows that one could nearly-cloak 
a delta point source inside the cloaked region. 

(ii) If one allows that the contents in the cloaked region could be arbitrary, 
then for a fixed near-cloak construction, there always exist cloaking-busting 
inclusions which defy any attempts of cloaking. These are similar to the res- 
onant inclusions observed in [9] for approximate acoustic cloaking. Following 
[9] , we employ a lossy layer with conducting medium outside the cloaked re- 
gion to overcome the resonance and re-achieve all the approximate cloaking 
results for passive media and active sources in (i). 

(iii) The performance of the approximate EM cloaking in asymptotically 
low and high frequency regimes. We show that it is impossible, with a fixed 
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near-cloaking scheme, to obtain cloaking uniformly in frequency, especially 
for the cloaking of active/radiating objects. Our observation is closely re- 
lated to the very recent study in [16], where frequency dependence for the 
approximate acoustic cloaking is considered. 

(iv) The limiting behaviors of solutions to regularized approximate cloak- 
ing problems, and their connections to finite energy solutions considered in 
[5] for singular ideal cloaking problems. 

Our study has been mainly restricted to spherical cloaking devices with 
uniform cloaked contents. We base our analysis on spherical wave functions 
expansions of EM wave fields. Nonetheless, we believe similar results would 
equally hold for general approximate EM cloaking study. 

In this paper, we focus entirely on transformation-optics-approach in con- 
structing cloaking devices. We refer to [4, 17, 23] for state-of-the-art surveys 
on the rapidly growing literature and many striking applications of trans- 
formation optics. But we would also like to mention in passing the other 
promising cloaking schemes including the one based on anomalous local- 
ized resonance [13], and another one based on special (object-dependent) 
coatings [1]. 

The rest of the paper is organized as follows. In Section 2, we present 
the basics on transformation optics in a rather general setting and apply 
them to the construction of EM cloaking devices. Sections 3-5 are devoted 
to the main results, respectively on, cloaking of passive media, cloaking of 
radiating objects and, cloaking-busting inclusions and retaining of cloaking 
by employing a lossy layer. The numerical experiments are given in Section 



2. Transformation optics and electromagnetic cloaking 

Let n be a bounded body in M 3 whose electric permittivity, conductiv- 
ity, and magnetic permeability are described by the M 3x3 - valued functions 
e, a and /x, respectively. Consider the time-harmonic electric field E and 
magnetic field H inside Q satisfying Maxwell's equations 



with uj > representing a frequency, J an internal current density. Let v 
be the exterior unit normal on the boundary dQ. By Ag a j we denote 
the linear mapping that takes the tangential component of E\qq to that of 
H\dQ, i.e., 



to recover the interior object, namely fJ.,£,cr and J, by knowing A^ j. It 
is pointed out that knowledge of the impedance map is equivalent to that 
of the corresponding scattering measurements (cf. [3]). We refer readers to 
[19] and [20] for uniqueness results of this inverse problem. Throughout the 
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rest of the paper, we shall denote by {Q,; e, /u, a, J} the object (EM medium 
and internal current) supported in Q. We would also use Aq to denote the 
impedance map in the free space; that is, it corresponds to the case with 
e = fj, = I,a = and J = in Q. In this context, invisibility cloaking can 
be generally introduced as follows. 

Definition 2.1. Let O and D be bounded domains in M 3 with D d $7. Q\D 
and D represent, respectively, the cloaking region and the cloaked region. 
{£l\D;e c ,[i c ,a c } is said to be an invisibility cloaking for the region D if 

K,»e,<re,J e = A o on dU for all u > 0, (2.3) 
where the extended object {O; e e , [i e , a e , J e } is given by 

;o . c „ „ n / e« Pc, 0c, 0} inOYD, 

{S2;e e ,// e ,<r e , J e } = < . (2.4) 

[ {£>;e a ,^ a ,cr a , J a \ mD, 

with {D;e a ,(i a ,o~ a ,Ja} being the target object (which could be arbitrary). 

According to Definition 2.1, the cloaking medium {f2\.D; e c , fJ-c, o~ c } makes 
the target object, namely the interior EM medium {D; e a , f-t a , o~ a } and the 
interior source/sink J, invisible to exterior boundary measurements. 

Next we present the transformation invariance of Maxwell's equations 
and transformation properties of EM material parameters, which shall form 
the key ingredients for our construction of invisibility cloaking devices. To 
that end, we first briefly discuss the well-posedness of the Maxwell equations 
(2.1). In the following, let Q be an open bounded domain in M 3 with smooth 
boundary. Assume that e, fj, and a are in L°°(J7) 3x3 , and they have the 
following properties: There are constants c m ,CM > such that for all x G £1 
and arbitrary ^6R 3 

Cm\H\ 2 <Z T e(x)Z<c M \Z\ 2 , c m \i\ 2 <e^x)i<c M \i\ 2 (2.5) 

and 

< i T ai < c M \C\ 2 . (2.6) 
We remark that the conditions (2.5) and (2.6) are physical conditions for 
regular EM media. We also assume that J G L 2 (fi) 3 . For the Maxwell 
equations (2.1), we seek solutions (E,H) E H(cwcl;{l) x H(curl; $7), where 

ff(curl; 0) = {u E L 2 (nf; V x u E L 2 (Q) 3 }. (2.7) 

We shall not give a complete review on the study of existence and uniqueness 
of solutions to (2.1) in the setting described above, and we refer to [14] for 
results related to our present study. It is noted that there is a well-defined 
continuous impedance map 

K,n,c,J ■ #~ 1/2 (Div; dQ) -> tf-^Div; dfl), (2.8) 

provided u avoids a discrete set of frequencies corresponding to resonances 
(cf. [14]). Here, 

H-^(Biv;dn) = {s E H~^(dtl) 3 ; s-u = a.e. on aOandDivs E H'^dQ)}, 
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with Div denoting the surface divergence on dQ. 

Lemma 2.2. Consider a transformation x = F(x) : Q — > Cl, which is 
assumed to be bi-Lipschitz and orientation-preserving. Let M = DF := 
(~^)ij=i ^ e the Jacobian matrix of F. Assume that (E,H) £ H(cnrl;Q) x 
-£f(curl; $7) are EM fields to (2.1), then for the pull-back fields given by 

E = {F~ X )*E := (M T y 1 EoF-\ H = (F~ X )*H := (M T )" 1 For 1 (2.9) 

and 

3= ^' j= Mm) mj ' f ~ 1 - i2M > 

we have (E,H) E i7(curl; f2) x 7/(curl;0) satisfying Maxwell's equations 

V x E = iujjlH, V x H = -iuj{e + i-)E + J in (2.11) 

LO 

where Vx denotes the curl in the x- coordinates, and £,fL,a are the push- 
forwards ofe,n,a via F, defined respectively by 

t = F,v:=^—M-v-M T oF-\ (2.13) 

* = F ^=dSM M '"- MT ° f '" (2 ' 14) 

Proof. The key ingredient for the proof of the lemma is the following trans- 
formation rule on curl operation (see, e.g. [14]) 

^ F = dSlM) M(Vx£)oF " ^ x " = J(M) M{VxH) ° F ~ K (2 ' 15) 

Using (2.15) along with (2.1), (2.9) and (2.13), we have 

V x E =— — — (V xfi)o F" 1 = , MUuifxH) o F~ x 
det(M) v ; det(M) v p ' 

=iu-^jj^M {iM T (M~ T H) o F- 1 = iuofiH. (2.16) 

Similarly, using (2.15), together with (2.1), (2.9), (2.10), (2.12) and (2.14), 
we have 

Vxfl = * (V x H) o F~ l = - —- (-iue r E + J) o F' 1 
det(M) y ' det(M) v 1 

= ~ iu - -— Me r E o F- 1 + - * MJ o F~ x 
det(M) det(M) 

= -iu£ r E + j, (2.17) 

where 

.o- , _ _ .a 

e r = £ + i— and e r = e + i—. 

LO LO 

The proof is completed. ■ 
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Corollary 1. Assume that F : f2 — > f2 is bi-Lipschitz and orientation- 
preserving with F\qq = Id. Using Green's identity, it is directly verified 
that 

v x E = v x E, v x H = v x H on dU, (2.18) 
which together with Lemma 2.2 yields 

K,!J,,a,J = ^e^^^aXF- 1 )*.}- (2.19) 

Lemma 2.2 and Corollary 1 summarize the basics of transformation op- 
tics in a rather general setting, which we shall make essential use of in the 
present paper. In the rest of this section, we give a short discussion on 
the singular ideal cloaking device construction considered in [5] and [18] 
using transformation optics, and introduce the notion of approximate cloak- 
ing from a regularization viewpoint. In the sequel, let B r denote the ball 
centered at the origin with radius r. Let M\ = B2, M 2 = B\ and M be 
the disjoint union M = M 1 U M 2 . Also, let N ± = B 2 \Bi, N 2 = B\ and 
N = N\ U N 2 . Moreover, set S := dB\. Consider the map 

Fx : Mx\{0} -> iVi, Fa(y) = M + l\y\\ ^ , < \y\ < 2 (2.20) 

which blows up {0} to N 2 while keeps the boundary dM\ fixed. In [5] 
and [18], the authors consider the lossless setting, i.e., one always assume 
that <t = 0. In the cloaking region Ni, the EM material parameters of the 
corresponding cloaking medium are given by 

(DF^HDFif 



«*)=*(*> = <*».'■ det(DFl) 



, x £ N\. (2.21) 

In the cloaked region N 2 = B\ , we consider cloaking an arbitrary but regular 
EM medium (eo,/io), he., 

jl(x) = n {x), e(x) = e {x) x £ N 2 , (2.22) 

which can be viewed as the push-forwards of (/xo, £0) i n M 2 by F 2 = Id. We 
denote the transformation by 

F = (F U F 2 ) : (Mi\{0},M 2 ) -> (iVi,iV 2 ). (2.23) 

By (2.21) together with straightforward calculations, we have in the standard 
spherical coordinates x 1— >• (r cos 4> cos 9, r sin (j> cos 8, r sin 8) that 

(r - l) 2 

/2 = g = 2 v ^ 2 y e r + 2e e , 1 < r < 2, (2.24) 

where e r and eg are respectively, the unit projections along radial and an- 
gular directions, i.e., 

j ~ ~T ~ "T " x 

It is readily seen that as one approaches the cloaking interface S the cloaking 
medium becomes singular, since i and jl no longer satisfy the condition (2.5). 
Finite energy solutions to the singular Maxwell's equations underlying the 
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cloaking are investigated in [5]. It is shown that {0} is a removable singular 
point. Specifically, let (E,H) be the EM fields corresponding to {N;i, /}}, 
then (E + ,H + ) = (Fi)*(E, H) are EM fields in free space on M\, which 
implies by Corollary 1 that 

AW _ AW 

A e,/i — A • 

On the other hand, (E~,H~) = (F2)*(E,H) satisfy the Maxwell equations 

V x E~ = iojp H~ , V x H~ = -iloeoE^ on M 2 , , 

v x B~ = 0, v x if- = on <9M 2 . ' 

Generically, one would have E~ = H~ = for (2.25) due to the homoge- 
neous 'hidden' PEC and PMC boundary conditions in (2.25) on 5M 2 . Also, 
due to such 'hidden' boundary conditions, it is claimed in [5] that one cannot 
perfectly cloak a generic internal current in the cloaked region B\. 

As can be seen from (2.24) the cloaking medium for the ideal cloaking is 
singular, which poses challenges to both mathematical analysis and physical 
realization. In order to construct practical nonsingular cloaking devices, it is 
natural to incorporate regularization by considering approximate cloaking, 
which we shall investigate in the subsequent sections. We conclude this 
section by introducing the notion of approximate EM cloaking. 

Definition 2.3. Let Q and D be bounded domains in M 3 with D <e Q, 
representing respectively the cloaking region and the cloaked region. Let 
p > denote a parameter and e(p) be a positive function such that 

e{p) ->■ as p ->■ 0+. 

{Q\D; £c, Pc, crc} is said to be an approximate invisibility cloaking for the 
region D if 

\\^M,J e - A o\\ = e(p) asp^0+, (2.26) 
where the extended object {Q; e§, o~§, J e } is defined similarly to (2.4) by 
replacing £ c ,p c , cr c with e p c ,p p c ,a p c . 

According to (2.26), with the cloaking device {Q,\D; e p , p p , we shall 
have the 'near-invisibility' cloaking effect. In order for the invisibility cloak- 
ing and approximate invisibility cloaking in Definitions 2.1 and 2.3 make the 
right sense, throughout the rest of the paper, we always assume that there 
is a well-defined impedance map Aq in the free space; namely, it is assumed 
that there is no resonance occurring in the free space. 

3. Nonsingular approximate cloaking of passive medium 

In this section, we consider the approximate EM cloaking for a relatively 
simpler case by assuming that all the EM media concerned are lossless, i.e. 
(7 = 0, and also there is no source/sink present, i.e. J = 0. 

For approximate acoustic cloaking by regularization, Kohn et al., in [9], 
proposed blowing up a small ball B p to B\ using a nonsingular transfor- 
mation F p which degenerates to the singular transformation F in (2.23) as 
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p — > 0, while Greenleaf et al., in [6], proposed blowing up B p to Br with 
R > 1 by the original singular transformation F. For the present study on 
approximate EM cloaking, we shall focus on the 'blow-up-a-small-ball-to-i?i' 
scheme and evaluate its performance. However, it is remarked that the other 
regularization scheme has been verified to yield the same performances for 
approximate EM cloaking. 

3.1. Construction of approximate EM cloaking. Let < p < 1 denote 
a regularizer and 



2-p ' 2-p 
Consider the nonsingular transformation from B<i to B<i defined by 



(3.1) 



. F^\y) = (a + b\y\)^ p < \y\ < 2, 
FJy) = { P m K ' V m, \v\ ^ 11 (3.2) 

' s 4Hv) = i \v\<p. 



Our approximate cloaking device is obtained by the push-forward of a homo- 
geneous medium in B2\B p by Fp 1 ^. Suppose we hide a regular but arbitrary 
uniform EM medium (eo,/io) in the cloaked region B\. Then the corre- 
sponding EM material parameter in B2 is 



(e p (x),p p (x)) 



((Fi%I,(F^%I) K|x|<2, 
(eq,Mo) \x\ < 1 



which are obviously nonsingular. The EM fields (E p ,H p ) € -ff (curl; .B2) x 
.ff (curl; #2) corresponding to {B2',e p , p p } satisfy Maxwell's equations 



V x E p = iujp p (x)Hp, V x Hp = —iue p (x)E p in B2, 
vxE p \ dB2 = f eH- 1 / 2 (Div;dB 2 ). 

By Lemma 2.2, the pull-back EM fields 

(E p ,H p ) = ((Fp)*Ep,(Fp)*H p ) G H(curl;S 2 ) x ff(curl;B 2 ) 

satisfy Maxwell's equations 

V x E p = iup p (y)H p , V x H p = -iu;e p (y)E p , in B 2 \B p , 
uxE p \ dB2 = f e H- l / 2 (Dw;dB 2 ), 

where 

( t\ ( w J ( J ' J ) p< |y| <2, 



(3.3) 



(3.4) 



By Corollary 1, we see that 



A w = A- - . 



Hence, the estimate of A~" ^ for the approximate EM cloaking is the same 
to that of A? „ . 
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3.2. Convergence and hidden boundary conditions. Henceforth, the 
following notations for EM fields shall be adopted 

E p := (E+,E;), H p := (H+,H~) for x G (B 2 \^, #i) 

and 

£ p := (£+ E~), H p := (H+, H~) for y G (B 2 \S^, 

We also use E 1 , i7 to represent the finite-energy EM fields considered in [5] 
for singular ideal cloaking which we discussed earlier in Section 2. (3.3) and 
(3.4) can be reformulated as the following transmission problems 



f V x E+ = iojfi p (x)H+, V x H+ = -iui p (x)E+ in B 2 \B U 
V x E~ = iunoH~, V x H~ = —iuiE§E~ in Bi, 
v x £+| E+ = ux E-\ s -, v x F+| E + = v x H-\x-, 
v x E+\ dB2 = f. 



(3.5) 



(3.6) 



and 

' V x E+ = iuH+, V x H+ = -iuE+ in B 2 \B~ P , 
V x E~ = iLjfi p (y)Hp, V x Hp = -iue p (y)Ep in B p , 

u x E p Ie+ = ^ x £p Ie- > " x ^Ie+ = " x ^ Is- » 
z/ x £?+| 9B2 = /• 

where S p := <9-B p . 

Our arguments rely heavily on expanding the EM fields into series of 
spherical wave functions. To that end, we introduce for n G Z + and m G Z, 

AQ(x) := V x {xj n (C|x|)r™0r)}, JV£ f (s) := V x {^(%|*|)>; m (z)}, 

where £ £ C is a complex number and x = x/|x| for i£l 3 , Here, Y™(x) are 
spherical harmonics and, hn\z) := j n (z) + iy n {z) with j n (z) and y n (z), for 
z G C, being the spherical Bessel functions of the first and second kind, re- 
spectively. The most important property of such functions for our argument 
are their asymptotical behavior with respect to small variables: 



M z ) 



0(\z\ n ), h n (z) = 0{\z\- n - 1 ), for |*| «1. (3.7) 



We refer to [3] and [14] for more properties of the functions introduced here. 

The second Maxwell's equations in (3.5) and the first of (3.6) would give 
rise to waves for x G B\ 



n=l m=—n /o c 

^ oo n \o.C 

^ = ^o 1/2 E E ^ 2 /3™M n % + a™VxM„%, 

n=l m=—n 
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and for y G B-2\B p 



oo n 



m 



E P = E E C n N n,u + CV X N™„ + j™M™ >0J + CV x M„™ , 

n=l m=—n 
1 oo n 

H p = ~ E E "'OC + CV x 2V£, + , 2 ^M™ + 7™V x M n " 

n=l m=—n 

(3.9) 

where /c = (pqEq) 1 / 2 . 

The following lemma characterizes the asymptotic behaviors of the coef- 
ficients in the spherical expansions (3.8) and (3.9) as p goes to zero. 

Lemma 3.1. Assume oj is not an eigenvalue of (3.5), namely, the corre- 
sponding homogeneous equations have only trivial solutions. Let (E p , H p ) 
be the unique solutions to (3.5), whereas (E p ,H p ) = ((F p )*E p , (F p )*H p ) be 
the unique solutions to (3.6). (E~,H~) and (E^,H^) are given by (3.8) 
and (3.9), respectively, whose coefficients are uniquely determined by the 
boundary data f . As p — > + , we have 

jZ = 0(1), Vn=0(l); c™ = 0(p 2n+1 ), d™ = 0{p 2n+1 ), (3.10) 

and 

a™ = 0(p n+1 ), (3™ = 0(p n+1 ). (3.11) 
Proof. We need to introduce the vector spherical harmonics 



CC := , = Grad F™, V'" : /' x !/> 
^n(n + 1) 



where Grad denotes the surface gradient. Define 

Unit) := h$(t) + thW'(t), J n (t) := j n (t) + tj' n (t). 

For t < 1, one can verify J n (t) = 0{t n ) and H n (t) = 0{t~ n - 1 ). Then on a 
sphere 8Br, we have for < R < 1, 

oo n 

uxE;\ dBR =e- 1/2 J2 E V^+l)(«n 3n(kuR)U™ 

n=l m=—n 

oo n 

ikuj 



1 oo n 

vx H-\ dBR = —p a l l 2 Y, E Vnin+^^kWjnikujR)^ 



n=l m=—n 

„m 1 



+a™j l J n (ku 1 R)V™ , 



(3.12) 
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whereas for p < R < 2, 



00 n 

uxE+\ 9Br = Y, E + 1) (<%hW(uR)U? + C-H n (uR)VZ' 

n=l m=—n 

+ 7 ™j n (u;R)U™ + ^MuR)V^ , 

1 oo n 

!/ x h+\ 9Br = ^ E E V^+iJPc^M)^' 



n=l m=—n 



(3-13) 

Expanding the boundary value on (9i?2 in terms of the vector spherical 
harmonics, we have 



oo n 



/ = EE (3.i4) 

n=l m=—n 

the boundary condition ^ x E£\qb 2 = f implies 
(R-l) 



c%h£\2u) + 1 ™j n (2u 1 ) = fW, 



Since E p = (F p 1 )*£'p, the transmission condition on the electric field in 
(3.5) reads 



v x E+\x+ = p(y x E+\z+) = vx E p 



By (3.12) and (3.13), we have 
(R-2) 



pc%h { n\up) + pj™j n (u;p) = e 1/2 a™j n (kuj), 
d™H n (iop) + «M = e~ 1/2 p™J n {kuj). 



Similarly, the transmission condition on the magnetic field implies 

(R 3) < fcc n ^n(wp) + kj™J n (up) = pQ 1/2 a™J n (ku), 
\ pC^ ) M+Wn^M=/Xo 1/2 Win(M- 

By (R-2) and (R-3), we have 

„m j. -,m „,m j. m jm + ^ra am + ^ra (n -i r\ 

c n - *l7n > «n - *27n ' d n ~ *3??n , P n ~ t ^n ) (3.15) 



12 APPROXIMATE ELECTROMAGNETIC CLOAKING 

where as p — > + , 

ti . = £ 1/2 kJ n (ujp)j n (ku) - fl~ 1/2 pj n (ujp)j n (kuj) = Qf p 2n+U 

P 1/2 ph ( i\ujp)J n (kuj) - e^ 1/2 k'H n (ujp)j n (kuj) 
^ . = kpj n (u)p)h ( h \ujp) - kpj n {up)rl n {ujp) = n+1 

p 1/2 pj n (kuj)hn\ujp) - e Q 1/2 kj n (kuj)ri n (up) lg \ 
^ . = Po 1/2 kJ n (up)j n (ku) - e 1/2 pj n (ujp)J n (kuj) = r p 2n+u 

e 1/2 ph£\ujp)J n (ku}) - p 1/2 kn n {ujp)jri{kuj) 

t4 . = pj n (ujp)h ( n\ujp) - pj n (Ljp)Un(uJp) = n+1) 

Eq 1 ^ 2 pj n (ku)hn\up) - Pq 1 ^ 2 kj n (kuj)'Hn(up) 

By (R-l), we have 

f W 9f (2) 

7 m = = Oil) n m = = O(l) 

ln t^^+uiu) U ' Vn t 3 n n (2u) + j n (2u) 

(3.17) 

By (3.15), these further imply 

a™ = 0(p n+l ), (3™ = 0(p n+l ), c™ = 0(p 2n+1 ), d™ = 0{ P 2n+1 ). 



We are in a position to evaluate the approximate EM cloaking. Our 
observations are summarized in the following. 

Proposition 3.2. For the approximate EM cloaking, if u is not an eigen- 
value of (3.5), we have 

l|AtAp- A oll=CV) asp^0+, (3.18) 

where \\ ■ \\ denotes the operator norm of the impedance map. 

Proof. We write the EM fields (E, H) propagating in the free space as 



oo n 



n = lr £T n n (3.19) 



IU! 

n=l m=—n 



Consider the boundary condition v x E\qb 2 = f satisfied by the (E, H) fields 
with / given by (3.14). By straightforward calculations, we have 

f (l) 9f (2) 
m _ Jnm , m _ ^Jnm 

n j n (2coY n J n (2u>Y 
Hence, the tangential magnetic field on the boundary is given by 

oo n , v 

iw { \ ILO 2 J 

n=l m=—n x ' 

(3.20) 
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Compared to ux H+\qb 2 from (3.13), one observes that c™, d™, 7™ — a™ and 
77™ — 6™ approach zero of order 0(p 2n+1 ), which in turn implies (3.18). ■ 

Proposition 3.2 shows that the approximate cloaking scheme constructed 
in Section 3.1 actually gives the near- invisibility cloaking effect. Next, we 
consider the limiting state of the approximate cloaking, showing that it 
converges to the singular ideal cloaking. 

Proposition 3.3. For the approximate EM cloaking, if u is not an eigen- 
value of (3.5), we have 

E p ^tE and H p -> H as p -> 0+. (3.21) 
Proof. We first show 

(E+,H+) = (FM)*(E+,H+)^(E,H) asp^0+. (3.22) 

It is easily verified that on any compact subset of B2 away from the origin, 
one has that {E+,H+) converges to (E,H) at the rate 0(p 3 ). Indeed, we 
shall show 

W E P ~ E Wl?(b 2 \W p ) + W H p ~ H Wl?(B2\W p ) = °(p 3/2 )> 
which implies (3.22). To that end, we note the following identities 

{ M^(x) = -yjn(n + l)j n {u\x\)V™(x), 
N^(x) = -y/n( n+l)hZ( u>\x\)V™(x), 

V x M%,(x) = ^ n( j" + l) J n {u\x\)U%(x) + n{n . + . l) jn{u\x\)Y?{x)x, 

\x\ \x\ 

V x N%Jx) = ^ x + 1) Hn(u\z\)U?W + ^^hWHx\)Y^(x)x. 

(3.23) 

By (3.9) and (3.19), we have 

00 n 



n=l m=—n 



+ ^ [« - + d%H n Wx\)]U™(x) 

+ n( ^[(C " K)3n{"\A) + d™hW{u\x\))Y™{x)x. 
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This implies as p — > + 
/ \Et-E\ 2 dx 



P 

'B 2 \B P 



oo n „2 

E E / TiCn + lJKTT-OinW + ^MlVdr 

n=l m=-n J P 



\ f 2 n(n + 1)| « - + d^Mfdr 

2 n 2 (n + 1) 2 |(^ - C)jnM + C^Wf^ 



E E °(/> 2n+1 ) = °(^ 3 



n=l m=—n 

by the convergence orders of the coefficients. Similarly, we have 



L 



p 

B 2 \B P 



H+ -H\ 2 dx = 0{p 3 ). 



On the other hand, it is observed from (3.12) that 

oo n 

E-\ s -=e 1/2 J2 E Vn(n + l)(a^j n (ku)U, 



v x 

n=l m=—n 

+f3™J n (kLu)V r 



III 
II 



m 

n I > 



oo 



^xf-| e - = ^- 1/2 £ E \A^+i)(^Vj„(M^n m 

+cC7„(fcw)lS 



n=l m=— n 



(3.24) 

are both of 0(p 2 ) as p — > + . Therefore the homogeneous PEC and PMC 
conditions naturally appears on the interior cloaking interface S _ . This 
is consistent with the interior 'hidden' boundary conditions discovered in 
[5] for singular ideal cloaking (see also (2.25)), which together with (3.22) 
implies (3.21). ■ 

3.3. Cloak-busting inclusions and frequency dependence. In our ear- 
lier discussion, we achieved near-invisibility under the condition that there 
are no resonances occurring. That is, u is not an eigenvalue to (3.3), or 
equivalently, to (3.4). In fact, if uj is an eigenvalue to (3.4), the small in- 
clusion (e p ,p p ) in the free space could have a large effect on the boundary 
measurement. In this resonance case, one may not even has a well-defined 
boundary operator A^ p ^ . The failure of the near- invisible cloaking due to 
such "cloak-busting" inclusion is also observed in the study of approximate 
acoustic cloaking in [9] . In the following, we shall show a similar result that 
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for a fixed approximate EM cloaking scheme, there always exists certain in- 
terior content (eo, no) such that resonance occurs at certain frequency oj. We 
shall be looking for such triples (oj, eo, po), or equivalently (oo, k, po), depen- 
dent on p, such that (3.5) is ill-posed. It corresponds to choices of (oo, k, po) 
such that the coefficient matrices of systems (R-l), (R-2) and (R-3) are 
singular. 

We consider two decoupled systems (R-l-l)-(R-2-l)-(R-3-l) and (R-l-2)- 
(R-2-2)-(R-3-2) corresponding respectively, to the variables {c™,a™,7™} 
and {d™, (3™, r]™}. The coefficient matrices are denoted as A n and B n in 
the following. By elementary linear algebra manipulations, the augmented 
matrix for A n for the first system reduces to 

/ h ( n 1) (2oj) j n (2u) f [ ni \ 



-e 



V 2 • (i \ ■ ( \ pj n (2ui)h£\up) j'jptfVp) 
Jn[KUJ) pj n \0Jp) (i), : (i) : 



V A n (3,3) A n (3,4) / 

where 

1/2 

4,(3,3) =^ty^ \p 1 /2 j n (koj)[J n (ojp)h ( n 1 \2oo) -H n (oop)j n (2oo)} 

h^ ] (2uj)j n (koo) 1 

-pp 1/2 J n (kuj)[j ri (u J p)h^(2oj) - h^(ojp)j n (2oo)]j , 

and 

in(M) = £ v 2f ™ \pp- 1/2 hW(up)J n (ku>) - pl /2 H n (Lop)j n (ku;)} . 

For det(A n ) = 0, one can choose (oo, k, po) satisfying 

jn(kuj) j n (ujp)h ( h\2uj) - h ( h\ujp)j n (2uj) 

^~n~n — \~P m ' \6.lv) 

MM j n (uj P )h { n\M - n n (ojp)j n (2uj) 

It is easily verified that with this choice, if fnm ^ 0, then 4,(3,4) ^ 0, 
there exists no solution of (c™,a™,7™). The boundary value problem is 
ill-posed and one does not have a well-defined boundary impedance map. In 
like manner, one can find (oo, k, eo) such that det(£? n ) = 0. 

Next we consider the performances of the approximate cloaking scheme 
in extreme frequency regimes. That is, we let p and (eo,po) be fixed, and 
evaluate the approximate cloaking effects as oj approaches zero or infinity, 
corresponding the low and high frequency regimes. First, we see that 

oo n 

" X H+\ BBa -UX H\ dB2 =J2J2 9nlUn + S^C. 

n=l m=—n 
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where 

_ -2iy/n(n + l)uj 2 t 3 b™W n (2uj) 

t 3 n n (2uj) + J n {2u) 
_ -2i v / n(n + l)u 2 b™ [e J n (ujp)j n (kuj) - pj n (ujp)J n {kuj)} W n {2uj) 

el ,2 det{B n ) 

9 n \ : =2^V»(n + l) (a^J n (2u) - J™Jn(2u>) - C% n (2uj)) 



J^n{n + l)t l a^W n {2uj) 

t 1 h < n ) (2uj)+j n {2u) 
_iy/n(n + ljgg [p J n (ujp)j n (kuj) - pj n {ujp)J n {kuj)} W n (2u) 

1/2 

fi ' det(A n ) 

(3.26) 

with W n (t) := j n (t)h n l) ' (t) - h£\t)j' n (t). 

We shall address the frequency dependence issue by assuming that the 
inputs are given by the EM plane waves of the form (6.1). The corresponding 

coefficients f nm and fnm are given by (6.2), while a™ and 6™ are given by 
(6.3). It is readily seen 

<C = 0(1), C = oiu- 1 ). 

For the low frequency regime with uj -C 1, by (3.26), it is straightforward to 
show 

(1) , n 2n+l (2) , ,n-1 2n+l 

9nrn ~ w P ) ffnm ~ w P 

which implies a satisfactory approximate cloaking. Whereas for the high 
frequency regime with w > 1, we exclude the influence of resonances from 
our study by considering the case that |det(-A n )| and |det(-B n )| are bounded 
from below by a positive function C nm (uj, p), where the transmission problem 
(3.5) is well-posed. Then we consider two separate cases: 

• When 1 < uj -C p^ 1 , i.e., ujp <C 1, since j n (t),h n . (t) oscillate be- 
tween — i _1 and t , J n {t) ,7i n (t) oscillate between —1 and 1, and 
W n (t) ~ t~ 2 as t increases, we have 

(1) (ujpyuj- 2 (2) (ujpruj- 3 
l9nml ~ c nm (u, P y l9nml ~ c nm (u, P y { - ' 

where one can show that C nrn (uj, p) < uj~ n ~ ?J p _n_1 . Here and in the 
following, for two expressions 1Z\ and IZ2, by "7£i < 7^2" we mean 
"Hi < cR-2 1 with a constant c. 

• For even higher frequency uj 3> 1 such that ujp > 1, we calculate 

| (1)|< — | ( 2 )|< — d 28) 

l9nml ~ c nm (uj, P y l9nml ~ c nm (uj, P ) { - ' 



APPROXIMATE ELECTROMAGNETIC CLOAKING 



17 



where C nm < w 2 . 

By (3.27) and (3.28), we cannot conclude whether or not one can achieve 
the near-invisibility. However, we conducted extensive numerical experi- 
ments to see that one has the near-invisibility in these two cases. These 
suggest that for the cloaking of passive media, excluding the resonance fre- 
quencies, one could achieve the near-invisibility for the approximate cloaking 
scheme of every frequency. This is sharply different from the case of cloaking 
active/radiating objects which we shall consider in the next section. 

So far, we have been concerned with the cloaking device where the cloaked 
region is B\ and the cloaking medium occupies B 2 \B\, which we obtain by 
using the transformation (3.1)-(3.2). It is remarked here that for arbitrary 
< R\ < R 2 < 00, one can construct an approximate cloaking device whose 
cloaked region is Br x and the cloaking layer is Br 2 \B\ by implementing the 
following transformation 



x ■= G p (y) 

where 



G { p\y) = (a + b\y\)ft p<\y\<R 2 , 
G?Hv) = * \V\<P, 



R\ — p _ R 2 — Ri 

--n-2) 



R 2 - P R2- P 

It is readily seen that all our earlier results hold for such construction. The 
remark applies equally to all our subsequent study. 



4. Approximate cloaking with an internal electric current at 

ORIGIN 

In this section, we consider the approximate EM cloaking scheme con- 
structed in Section 3.1 in the case that we have an internal electric current 
present in the cloaked region supported at the origin. The corresponding 
EM fields verify 

V x E p = iujppHp, V x H p = -iujEpEp + J, in B 2 ^ 
v x Ep\ dB2 = /, 

where J has the form 

J= E (W(*))v aj (4.2) 

\a\<K 

with 5o denoting the Dirac delta function at origin and v a E C 3 . The 
pull-back EM fields satisfy 

V x E p = iujppHp, V x Hp = —iujEpEp + J, in B 2 , ,^ ^ 
v x Ep\ dB2 = /, 

where J = (F p 2 )* J. 
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The point electric current J would give rise to a radiating field 

K n 

E J = H E Pn N n >ku + CV X N™ huy (4.4) 
n=l m=—n 

Hence for x £ B\ 



V? = £ E E «™ M n% + /CV x Jl^ + p™N™ kuJ + q™V x A^, 
n=l m=— n 

where and q™ equal zero when n > K. Whereas E+ and are as in 
(3.9). 

Lemma 4.1. Assume u is not an eigenvalue to (4.1). Let (E p ,H p ) be the 
EM fields to (4.1) with J given by (4.2), and (E p ,H p ) = ((F p )*E p , (F p )*H p ) 
be the EM fields to (4.3). Given E~ as in (4.5) and E+ as in (3.9), we have 
as p — > + , 

7™ = 0(1), C = <?(1); C = 0( P n+1 ), d™ = 0( P n+1 ), (4.6) 

and 

a? = 0(1), ff? = 0(l). (4-7) 

Proof. The boundary condition on <9-E?2 implies (R-l). From the standard 
transmission conditions, we have 

(X%h£\u, P ) + PlniniuP) = ^ 1/2 (*n(M +p^ 1) (fco;)). 

d%H n (u;p) + vZMup) = e- 1/2 (p™J n (kcj) + q™H n (kuj)), 

kc™U n (up) + hyZJ n (up) = p^ ll2 {a^J n (ku)+p^n n (ku)), 
pd™h£\ujp) + prtfjniup) = p 1/2 (kf3™j n (koj) + kq%h£\ku)). 

Solving (R-2') and (R-3'), we obtain 

„m + -,m i j./ _.m _,m j. ,,m _j_ j./ „m 

c n - *l7n + *lPn > a n ~ *27„ + *2Pn . f 4 s>\ 

«n — *3 r /n + t 3^n ' Pn ~ t i r ln + t i ( ln , 

where (i = 1, 2, 3, 4) are given by (3.16) and t[ (i = 1, 2, 3, 4) are given by 

f , = h { n ] (ku)J n (ku) -H n (ku)j n (ku;) = 0{p n+l ) 

P 1/2 phn\ujp)Jn(kuj) - e 1/2 k'H n (ujp)j n (kuj) 
t , _ EQ 1/2 kh i h\kuj)'H n (ujp) - p 1/2 p'H n {kuj)h ( n\ujp) _ 

p 1/2 ph ( h\ujp)J n (kuj) - e Q 1/2 k'H n {ujp)j n (kuj) ^ ^ 




J n (ku)h^\kuj) -n n (kto)j n (kuj) _ Q, n +l\ 

e 1/2 phn\ujp)J n (kuj) - Pq 1/2 kU n {up)3n{ku) 



t , P 1/2 kh < h\ku))H n (up) - £ 1/2 pH n (kuj)h < n\ujp) _ 
e 1/2 ph ( n \ujp)J n (kuj) - Pq 1/2 k% n (uj p)j n (kuj) 
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Plugging into (R-l), we obtain 

m _ fnm ~ P™t'\hn _ (2cd) _ , . m _ 2fnm — t' 3 q™'H n (2u)) _ . 

ln t^w+MM t 3 n n (2u) + M2u) uu ' 

(4.10) 

which together with (4.8) imply (4.6) and (4.7). ■ 

Next, we evaluate the performances of the approximate EM cloaking. 

Proposition 4.2. For the approximate EM cloaking with an internal point 
current (4-2) present in the cloaked region, if uj is not an eigenvalue to (4.1), 
we have 

\\k" Sp!i . pJ -k%\\=Q{p 2 ) asp^Q+, (4.11) 

where || • || denotes the operator norm of the impedance map. 

Proof. On the boundary dP>2, using the expression (3.13) for v x H+\qb 2 
and (3.20) for v x H\qb 2 , together with the asymptotic estimates of the 
corresponding coefficients in Lemma 4.1, we have (4.11) by straightforward 
comparisons, since the coefficients c™, d%, 7™ — a™ and n™ — b™ converge 
to zero of order 0(p n+1 ). ■ 

By Proposition 4.2, we see that one still achieves near-invisibility cloaking 
even though there is a source/sink present in the cloaked region. That is, 
the approximate cloaking makes both the passive medium and the active 
point source/sink nearly-invisible. However, we have one order reduction of 
the convergence rate. This is due to the extra terms 

-p™t' x h£X2u)) -q%t' 3 H n (2uj) 



tih { n\2L0) - j n (2u) ' *3^n(2w) - J„(2w) 



, «, ~ p 



n+1 



in 7™ — a™, r\™ — b™, c™ and d™ respectively, compared to the case without 
the source/sink. 

Next, we consider the limiting status of the approximate cloaking in this 
case when a point source/sink is present. We have 

Proposition 4.3. Assume uj is not an eigenvalue to (4.1). Let (E p ,H p ) be 
the EM fields satisfying (4.1) and (E p ,H p ) = ((F p )*E p , (F p )*H p ) be the EM 
fields satisfying (4.3). Then we have as p — >• + , 

(E+,H+)^(E,H) (4.12) 

with (E,H) being the EM fields on B2 in the free space. Also 

(E^H;)^(E-,H-), (4.13) 

where (E~,H~) satisfy the Maxwell equations 

V x E~ = iu)n H~, V x H~ = -iuje E~ + J in B x (4.14) 

with 

z,xir| s -^0 and uxH~\ E -^0. (4.15) 



20 APPROXIMATE ELECTROMAGNETIC CLOAKING 

Proof. By a similar argument to the first part of the proof of Proposition 
3.3, one can show that on any compact subset of B2 away from the origin, 
(E+,H+) (E,H) at the rate 0(p 2 ), and on B 2 \Bp, 

W E P ~ E \\lHb 2 \b- p ) + \\ H p ~ H WmB 2 \B p ) = 0{p 1/2 ) as p -+ 0+. 

This proves (4.12). Next, we shall show (4.15) which in turn implies (4.13)- 
(4.14). On the interior cloaking interface S - , the Cauchy data are given 
by 

00 n 

uxE-\^-=e- 1/2 J2 E v / ^+i)((<in(M+rf(M)C 

n=l m=—n 

+(f3™J n (kco) + q™H n {kuj))V™) , 

— 1/2 00 n 

» x h-\ s - = f^- E v / ^+i)(« ^n(M +p™ ^(m)C 1 

n=l m=—n 

+&J*{ffl3n{ku) + (£l£\ku>))U™). 

(4.16) 

We observe that as p — >• + 

/3™Jn(M + C?*n(M = U^J n (ku) + (t' 4 J n (koj) + U n {kuj))q™ = 0(1), 

(4.17) 

where 

+' • ri. \ , ^ p$\up)\j n {ku))'H n (ku) - fe^ 1) (fcq;) > 7 n (fctj)] 

t 2 j n {ku)) + h\>{ku}) — — = (p), 

p Ji n {u}p)j n {ku) 

+' <T fu fu \ jn{ku)rL n {ku) - J n {ku))h£\ku) 

t^Jnykuj) + U n {ku) — r = 0(1). 

Jn{kuj) 

Similarly, we have 

c%J n (ku>) + P ™U n {kLo) = 0(1). (4.18) 

Plugging (4.17) and (4.18) into (4.16), we have (4.15). The proof is com- 
pleted. ■ 

By Proposition 4.3, we see that as p — > + , the near-cloak converges to 
the ideal-cloak. Moreover, in the limiting case, the EM fields in the cloaked 
region are trapped inside and the cloaked region is completely isolated. 

Finally, we consider the frequency dependence for the approximate cloak- 
ing of active/radiating objects. Again, we address our study by considering 
the inputs being EM plane waves as in Section 3.3. By straightforward cal- 
culations, the coefficients that characterize the difference of the boundary 
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measurements, i.e. v x H~^\qb 2 — v x H\qb 2 -, associated to terms U™ and 
V™ verify 

gW :=%Jn(n + T){b™j n {2uj) - r]™j n (2u>) - d^{2u:)) 



i 

-2iy/n(n + l)oj 2 W n (2uj) 

el /2 det(B n ) 
-el /2 q™kojW n (koj] 



b™ (e J n (ojp)j n (kuj) - pj n {ujp)J n (ku)) 



9 { ni :~Vn(n+l)(a™J n (2u) - ^J n (2oj) - c™U n {2u)) 
i^n{n + l)W n {2uj) 



1 II 

p ' det(A n ) 
f/ /2 p™ku;W n (kuj) 



C (PoJn(up)jn(kuj) ~ pj n {ujp)J n (ku)) 

(4.19) 



In the low frequency regime with uj <C 1, by (4.19) we have 

which implies that one cannot achieve near-invisibility when ui < p 2 / 3 . In 
the high frequency regime with w > 1, by excluding the resonances and 
using similar arguments to that in Section 3.3, one can show 

|5« I < |o (2) I < (4 20) 

where 

u -n-3 p -n-l up < 1, w > 1, 



C nm (w,p) 



< 



w 2 > 1. 



By (4.20), one cannot conclude whether or not the near-invisibility is achieved. 
However, in our numerical experiment given in Section 6.3, we have observed 
the failure of the approximate cloaking in the high frequency regime. There- 
fore, it can be concluded that for a fixed approximate cloaking scheme with 
a point source/sink (4.2) present in the cloaked region, in addition to reso- 
nances, the near-invisibility cannot be achieved uniformly in frequency. 



5. Approximate cloaking with a lossy layer 

In our earlier discussion of lossless approximate cloakings, we have seen 
the failure of the near-invisibility due to resonant inclusions. Following the 
spirit in [9] by introducing a damping mechanism to overcome resonances in 
approximate acoustic cloaking, we surround the cloaked region first by an 
isotropic conducting layer, then another anisotropic nonconducting layer as 
described as earlier. To be more specific, given a damping parameter r > 0, 
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our new regularized parameter in B 2 is given by 

f {{F 2p )J 1 {F 2p )J) K|x|<2, 
(fi T ,e T ) := ((F 2p )*I,(F 2p )*(l + iT)) \ < \x\ < 1, 



(p p (x),£ p (x)) = < 



, (Mo,£o) 



\x\ < 9 , 



which is the push-forward of 

' (1,1) 

(v P (y),£ P (y)) = < 



2p < \y\ < 2, 
{1,1 + it) p<\y\<2p, 
[ ((F^)^ ,(F^s ) \y\< P . 



(5.1) 
(5.2) 



by the transformation 
x := F 2p (y) 



(^ + 2(T^)M)ft 2p<H<2, 

i \v\< 2p. 



To assess the approximate cloaking in this setting, we consider the trans- 
mission problem 

' V x E\ = iuHi, V x Hi = -iuEi, in 2p < \y\ < 2, 

V x E 2 = iup T H 2 , V x H 2 = —icu£ T E 2 , in \ < \x\ < 1, 

V x E 3 = icjpoHz, V x H 3 = -iueoEs + J, in \x\ < \, 
v x Ei\ dB2 = /; 



" x E i\dB- = 2 P( V x ^i)Isb+ ' v x #2^- = 2p{u x fTi)| aB + ; 

(5.3) 

The problem is well-posed on B 2 since e T is complex. Actually, we have 

(p T ,e T ) = (2p, 2p(l + ir)). 

Set 

/c T := (^ T e r ) 1/2 = 0(p) as p -> 0+. 
We can write the spherical wave expansions of the electric fields as follows 



00 n 



n=l m=—n 

00 n 



n=l m=—n 
00 n 



^3 = £o V2 £ E < M n% + W V X + p^iV^ + <?™V X N^. 

n=l m=—n 

(5.4) 

Then we have 

Proposition 5.1. For any uj € R + , assume the EM field (E p ,H p ) satisfies 

j V x E p = iujp p H p , V x H p = —iioipEp + J, in i?2 
\ v x £ p |aB 2 = /, 
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where {fi p ,i p ) is the lossy medium given by (5.1). Then the pull-back field 
(E p , Hp) satisfies 

V x E p = iuj^pHp, V x Hp = —iujEpEp + J, in B2, 
v x E p \ dB2 = f 

with (fip,e p ) given by (5.2). Therefore, the fields 

(Ei, Hi) = (E p \b 2 \b^> h p\b 2 \W;)> 
(E 2 ,H 2 ) = (E p \ Bx ^ i2 ,Hp\ BA - g ^), 

(E3,H 3 ) = (Ep\ Bl/2 ,Hp\ Bl/2 ) 

satisfy the transmission problem (5.3). Moreover, 

(i) If J = 0, then (Ei,E 2 , E 3 ) is given by (5.4) with p™ = q™ = for all 
n and m, and 

7 ™ = 0(i), 77™ = 0(1), c = 0(p 2 «+ 1 ), c = 0(p 2n+1 ); 
C = 0(/> 2n+5/2 ), C = 0(p 2n+5/2 ), 7^ = 0(p 3/2 ), C = 0(p 3/2 ); 

< = 0(p" +1 ), /3™ = 0(p n+l ). 

(5.5) 

(ii) If J ^ Q is given by (4.2), i/ien (Ei,E 2 ,E 3 ) is given by (5.4) wrai/i 
PrT) 9™ 7^ / or some n an<i m, and 

7^ = 0(1), *C = 0(i), C = 0( P " +1 ), c = 0(p n+1 ); 

C™ = 0(p"+ 3 / 2 ), = 0(p n + 3 / 2 ), 7™ = 0(p- n+1 /2) fj% = 0(p~ n+1 /2) ; 

a- = 0(1), ^ = 0(1). 

(5.6) 

Proof. In the case that no source/sink is present (J = 0), the boundary 
condition and transmission conditions in (5.3) imply (R-l) and the following 
equations. 

e- 1 ' 2 + <Zhg>(krU,j\ = 2p [^jn{2up) + c™h£\2u;p)) , 

Er V2 (fj^MKcj) + dZHnikrU)) = VnM^P) + d™Hn(2up). 

Pr 1/2 (CU n (k T u) + i™J n {k T u)) = k T (CU n (2cup) + j™J n (2cjp)) , 
Pr 1/2 k T (d%h { n\k T 0j) + r,™j n (k T Uj)) = 2p (d^h^(2up) + r]™j n (2ujp) 




^0 1/2 «nin(^) = Sr 112 ^j n {ht)+C^\ k -^)) , 

-- 1 / 2 am rr /fej\ _ g.^ 1 / 2 / z.m rr ( k T uj \ , Jm-iy ( k T uj \ 
'0 Jn{ 2 ) — t-r I 'In 2 > ' n rL n\ 2 J 



(R-6) 



, , f Mo V2 *V< J n (f) = Pr ll2 k (c-^ n (Y) + 7»(^)) > 
1 " j \ /^ 1/2 Wj„(^) = Av/^ 1/2 (C^ 1) (Y) + &(Y))- 
Solving (R-6-1) and (R-7-1), we obtain 

g m = / ia m ^ m = 7 oQ ; m 
C„ — l\VL n , J n — l 2 U- n , 
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where as p — > + 

j _ £ ^ g _ (jn(~2~)t7n(~2~) ~ PtPq i7n(~2~) jn(~2~)) _ ^, n +3/2\ 



- C(p _n+1/2 ) 

Plugging the above quantities into (R-4-1) and (R-5-1), we further have 
(R-4-1) - na: + 2p £ y 2 / l < l 1 W)C = -2 /3e y 2 i n (2^ /) ) 7 -, 
(R-5-1) - r 2 a™ + p 1 J 2 k T n n (2ojp)C = -^ /2 fcr JnpuphZ, 

where 

ri = hh£\k T u) + Z 2 j n (^) = ©(p 1 / 2 ), 
r 2 = hH n (k T uj) + l 2 J n {k T Lo) = 0(p 1/2 ). 

Then 



c m = sTy m a m = soY™ 



where 



_ Mr J n (2up)ri - 2pj n (2ujp)r 2 _ 2n+1 
2ph} n ' (2ujp)r 2 - p T Un{2u)p)ri 

2pp T eJ 2 (-j n {2up)'H n (2ujp) + J n (2ojp)h i n ) (2copy 

S2 = ^ m L = 0(p n+u 

2phn'(2ujp)r 2 - p T 'H n {2ujp)ri 



ln= m . — = 0(1), (5-7) 



By (R-l) as /> — >■ + , we have 

r(l) 
nm 

s 1 h ( a\2co)+j n (2u;) 

which in turn implies 

C = 0{p 2n+1 ), a™ = 0(p n+1 ), 

c™ = <D(p 2n+5 / 2 ), ^ = 0(p i ' 2 ). 
Similar calculations suggests the other estimates in (5.5). 

Statement (ii) is derived from solving (R-l), (R-4), (R-5) and 

£0 (jnC-rX + t£\%)l%) = Zr 1/2 (&(Y) + C^ } (^)) , 

e 1/2 (J n (f )/3™ + W n (^)C) = ^ 1/2 (TO( Y) + C«n(¥)) > 
M o 1/2 fc r (<\7n(^) +p^7£ n )) = ^ V2 A: (c£« n ( ^) + 7^«( ¥)) . 
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Using Proposition 5.1, all our results in Sections 3 and 4 for the loss- 
less approximate EM cloaking can be shown to hold equally for the lossy 
approximate cloaking scheme (5.1). We remark briefly on this here. 

Remark 5.2. The estimates in (i) of Proposition 5.1 imply that, without 
an internal source/sink, the EM fields (-£3,^3) in the cloaked region £?x/2 
degenerates in order 0(p 2 ). Whereas for the EM fields (E2, H2) in the lossy 
layer Bi\Bi/ 2 , it is easily seen 

Then (5.5) and (5.4) imply that (E 2 ,H2) degenerate in order 0{p 2 ) as p 
decays. It follows that the vanishing Cauchy data appears on the inner 
surface dB^ . Moreover, the boundary operator on dB 2 of the approximate 
cloaking converges to that of the ideal cloaking in order 0(p 3 ). 

Remark 5.3. With an internal point source/sink of the form (4.2) present 
in the cloaked region, the asymptotic estimates of the coefficients for the 
corresponding EM fields are given in (ii). By straightforward verification, 
one can show near-invisibility for the lossy approximate cloaking similar to 
Proposition 4.2 in the lossless case. On the other hand, one can also show 
that both the EM fields (£2, #2) and (£3, #3) are 0(1), and hence they 
do not degenerate. Moreover, the Cauchy data on the inner surface dB^ 
does not vanish since by (R-4) and (R-5), the terms associated to V™ of 
v x E 2 and v x H2 are 0(1). These observations suggest that in the limiting 
case, the lossy approximate cloaking converges to the ideal cloaking, and 
the cloaked region is completed isolated with the EM fields trapped inside 
(see Proposition 4.3 for similar observations in the lossless case). 

For the frequency dependence of the performances of the lossy approx- 
imate cloakings, we also have completely similar results to those in the 
lossless case, which we would not repeat here ( see our discussion at the end 
of Sections 3 and 4). 

We conclude this section with two more interesting observations. In [9], 
for the approximate acoustic cloaking by employing a lossy layer, one needs 
to require that the damping parameter r ~ p~ 2 , which is not necessary for 
our present approximate EM cloaking. On the other hand, it is shown in [16] 
that if r is allowed to be p-dependent, one could achieve near- invisibility uni- 
formly in frequency. However, such result does not hold for the approximate 
EM cloaking. 

6. Numerical experiments 

In this section, we carry out some numerical experiments based on the 
discussions and calculations in Sections 3, 4 and 5. First we introduce an 
electric plane wave of the form 

E = e ~ iuJX - d P (6.1) 
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with d = (1, 6 d , 4> d ) € § 2 , P G C 3 and d ■ P = 0. In the free space, the EM 
fields (E, H) := (e~ iux - d P, -e~ iwxd d x P) satisfy Maxwell's equations 

VxE = iujH V x H = -IujE. 

The spherical wave functions expansions of the EM-fields (E, H) are given 
by 

( oo n 

£ = EE «;(x) + cvxM3 

n=l m=—n 
1 oo n 

^=-EE "'CM™ (*) + <VxM n ™ (x), 



n=l m=—n 



where 



and 



/n 
jn(2w)' 



2/A 

Jn(2w)' 



f (l) 

f (2) 
Jnm 



f nm (d,P) 

f£L(d,P) 



4tt 



n[n + l)i n 

47T 



M™(2d)-P, 



n(n + 



^-V x M%j2d) ■ P. 



(6.2) 



By (3.23), we have 

4-7T 



4tt 



On the boundary dP>2, one has 



U™{d) ■ P. (6.3) 



x x e -^' d p 



n(n + \)ui n 

8b 2 = Y, E v / ^Ti)(/^tC(^) + /^iC(^ 

n=l m=— n 



(6.4) 

Figure 1 demonstrates an electric field by taking the first 15 modes in 
the above expansion; that is, n is up to N = 15. Throughout all our 
computations, we shall make use of such truncation when a spherical wave 
function expansion is considered. 

6.1. Lossless approximate cloaking of passive media. Recall in Sec- 
tion 3 that the EM material parameters of our lossless cloaking device are 

((ifV,(^V) K|x|<2, 
(eoiMo) ~ arbitrary constant \x\ < 1. 



(i p (x),fLp(x)) 



Based on the calculations in Lemma 3.1, we depict the EM fields (E p ,H p ) 
propagating in {B2',i p , fi p } in Figure 2 with the following boundary condi- 
tion 

i x E p \qb 2 = x x E, (6.5) 
where E is the one demonstrated in Fig 1. It is remarked that the boundary 
input (6.5) will also be implemented in our subsequent numerical exper- 
iments, when a boundary condition is concerned. Next, we consider the 
convergence of the near-cloak to the ideal-cloak. To that end, for the EM 
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Figure 1. Real part of E±, namely the first component of 
E (sliced at x = 0, 1, 2), with the first 15 modes and to = 5, 
d = (l,7r/2,7r/2) gS 2 ,P = (1,0, 0) T . 




Figure 2. Real part of (E p )i (sliced at x = 0,1,2), with 
oj = 5, e Q = hq = 2, p = 1/6. 



fields (Ep,H p ), we compute the deviations of the boundary operators via 
the formula 

Er(p) := \\x x H a — x x H\\ i„. 
In our calculations, we shall make use of the following identity from [14], 



oo n 



H~2(DiV;dB2) 



Yl S Vn(n + l)\g ( ^ m \ 2 + 



n=l m=—n 



- v /n(n + 1) 



k (2) I 2 

|yn,ml i 
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given the vector spherical harmonic expansion of A 

oo n 

A = V V U m + V m 

n=l m=—n 

The convergence rate as p — > + is calculated as following 



Er(pi) pi 
r{p) := In / In—. 

Er{p 2 )l p 2 



Pi,P2 -> + . 



In Table 1, we list the computational results, which verify Proposition 3.2, 
i.e., the convergence order is 3. 



p 


0.1 


0.05 


0.01 


0.005 


0.002 


0.001 


Er(p) 


0.1810 


0.0139 


8.42e - 05 


1.02e-06 


6.42e - 07 


7.97e - 08 


r{p) 




3.703 


3.173 


3.044 


3.020 


3.009 



Table 1. Convergence rate of boundary operator for the 
lossless approximate cloaking with u = 5, £q = po = 2. 



6.2. Lossless cloaking of active/radiating objects. In this numerical 
experiment, we study the performance of our lossless approximate cloaking 
device when an internal point source/sink is present at origin, elaborating to 
the discussion in Section 4. We apply a delta source J = Y1\u\<k(^x ^o(^)) v a 
by introducing a radiating field 

K n 
n=l m=—n 

with known and , into the electric field E~ inside the virtue inclusion 
Bp. In Figure 3 and Table 2, we choose p^ 1 = p\ = p\ = 5, q^ 1 = q® = q\ = 
2 and = p™ = otherwise. From Figure 3, we see that one could still 
achieve near-invisibility, and the EM fields in the cloaked region is almost 
trapped inside. Table 2 verifies that the convergence order of the near-cloak 
is 2, which is consistent with Proposition 4.2. 



p 


0.1 


0.05 


0.01 


0.005 


0.002 


0.001 


Er{p) 


1.9787 


0.3509 


0.0114 


0.0028 


4.41e-04 


1.10e-04 


r(p) 




2.495 


2.129 


2.031 


2.013 


2.006 



Table 2. Convergence rate of boundary operator for the 
lossless approximate cloaking with a delta source, w = 5, 
£o = P-o = 2. 
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Figure 3. Real part of (E p )i for the cloaking problem 
(sliced at x = 0, 1, 2) with a delta source at the origin, oj = 5, 
£o = = 2, p = 1/12. 

6.3. Cloak-busting inclusions and frequency dependence. In Section 
3.3, we have shown the failure of lossless cloaking due to resonances. In 
Figures 4, for a fixed p, the first mode (n = 1) of boundary errors Er(p) 
are plotted vs frequency to, for both passive and active cloaking. We ob- 
serve blowups of the errors at resonant frequencies, where the determinants 
det(^4 n ) and det(B n ) (n = 1) vanish (see Figure 5 for those resonance fre- 
quencies). In fact, we have numerically shown that for every frequency co 



Lossless passive cloaking (p-0.01) 



Lossless active cloaking (p-0.01) 




Figure 4. Boundary error for mode n = 1. Left: lossless 
cloaking (no source). Right: lossless cloaking (with a source). 
p = 0.01, u G [1,3]. 



and p, there is a choice of 'cloaking-busting' inclusion in B\, e.g., a pair of 
parameters (£q,/xo) satisfying (3.25), such that the lossless construction is 
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_ 30 
< = 

2 20 




Figure 5. Frequency dependence of determinants of coef- 
ficients system (R-l), (R-2) and (R-3). p = 0.01, oo G [1,3]. 

resonant. In Figure 6, an example of such resonant inclusion at mode n = 1 
is plotted against p for a fixed frequency. One can see that as p — > + , 
the EM parameters of the inclusion become singular, namely, Eq — > oo and 
fiQ — > 0. As discussed in Section 3.3 and Section 5, Figure 7 demonstrates 



Cloak-busting inclusion at 01=14, n=l, k=l 
25 1 1 1 1 . 1 2.5 1 1 




0.02 0.04 0.06 0.08 0.1 0.02 0.04 0.06 0.08 0.1 

P P 



Figure 6. EM parameters (pq,£q) for a cloak-busting in- 
clusion at n = 1 cj = 14 with k = 1. 

that both the lossless (excluding the resonant frequencies) and lossy cloak- 
ing schemes work well in the low frequency regime, namely when uj <C 1, 
without any source/sink present in the cloaked region. In Sections 4 and 
5, when a point source/sink is present at the origin, we see that both the 
lossless and lossy cloaking schemes fail when oj < p 2 ^ 3 , as shown in Figure 
8. For higher frequencies, the behaviors of the cloaking schemes are not 
deterministic. Nonetheless, we show in Figure 9 that the lossless cloaking 
of active/radiating objects (excluding resonant frequencies) generates rela- 
tively large boundary error Er{p) when cj 3> 1. 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 I 



Figure 7. Approximate cloaking performance in low fre- 
quency regime uj 6 [0,1]. Left: boundary error (n = 1) for 
lossless cloaking (no source). Right: boundary error (n = 1) 
for lossy cloaking (no source), p = 0.01. 

Lossless active cloaking in low frequency regime (p-0.01) 
250 1 — I 1 1 1 1 



200 



150 

Q. 

100 




0.05 0.1 0.15 0.2 



Figure 8. Boundary error (n = 1) for cloaking with a 
source. Approximate cloaking compromises in low frequency 
regime, p = 0.01. 

6.4. Lossy approximate cloaking. According to our discussion in Section 
5, we employ a lossy layer right between the cloaking layer and the cloaked 
region. In Figure 10, we show how the EM-fields propagate in such a lossy 
construction of approximate cloaking. One can see that near-invisibility is 
achieved. In Table 3, the convergence order of the lossy near-cloak of passive 
media as p — > + is shown to be 3, which is consistent with Remark 5.2. It 
is recalled that for the lossy approximate cloaking, the EM parameters in 
B 2 are given by 

( {(F 2p )J,(F 2p )J) 1 < \x\<2, 

(fi p (x),i p (x)) = < (p T ,e T ) := ((F 2p )J,(F 2p )*(l +ir)I) \ < \x\ < 1, 
I (a*o,£o) M < |- 
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Lossless active cloaking in high frctjucncy regim (p-0.01) 




1000 1000.5 1001 1001.5 1002 1002.5 1003 1003.5 1004 1004.5 1005 



Figure 9. Approximate cloaking (with a source) perfor- 
mance in high frequency regime . Boundary error (n = 1) 
Er{p) > 2 for uj £ [1000, 1005]. p = 0.01. 



At last, we demonstrate the frequency dependence of our lossy approximate 




Figure 10. Real part of {E p )\ for the lossy approximate 
cloaking problem (sliced at x = 0,1,2), with ui = 5, eq = 
po = 2, p = 1/6. 



p 


0.1 


0.05 


0.01 


0.005 


0.002 


0.001 


Er{p) 


0.2733 


0.0455 


3.75e - 04 


4.69e - 05 


3.00e - 06 


3.75e - 07 


r(p) 




2.5867 


2.9818 


2.9998 


2.9998 


2.9998 



Table 3. Convergence rate of boundary operator for lossy 
approximate cloaking of passive medium, with frequency uj = 
5, £q = Mo = 2, damping parameter r = 3. 
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cloaking scheme in Figure 11 without a source/sink. Observe that the reso- 
nant frequencies disappear. However, we observe some frequencies at which 
the boundary error Er(p) is relatively large. We believe such frequencies 
are those very close to the poles or transmission eigenvalues in the complex 
plane of the boundary value problem. It is remarked that such phenome- 
non could also be observed in the lossless approximate cloaking. If there is 
a point source present at the origin, we would have the similar numerical 
result as the case considered in Figure 9 for the lossless cloaking. 

Lossy passive approximate cloaking (p-0.01) 



1 


I. . 


. 1 , 


1 1 J 




, J 


I 



123456789 10 



Figure 11. Boundary error (n = 1) of lossy approximate 
cloaking (no source), p = 0.01. 
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